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Abstract 

In this paper we show how to calculate explicitly the exponential of certain matrices, 
which are evolution operators governing the interaction of the four level system of atoms 
and the radiation, etc. We present a consistent method in terms of the magic matrix by 
Makhlin. 

As a closely related subject, we derive a closed form expression of the Baker-Campbell- 
Hausdorff formula for a class of matrices in SU(4), by use of the method developed by 
the present authors in |quant-ph / 06 10009 , 



1 Introduction 

The purpose of this paper is to develop a useful method to calculate the exponential of certain 
matrices explicitly. These matrices arise from differential equations governing the interaction 
of the four level system of atoms with the radiation whose image is illustrated in Figure 1. 



*E-mail address : fujii@yokohama-cu.ac.jp 
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Fig.l Atom with four energy levels and general action by laser fields 

For the general background, motivation and possible applications of the present work, see the 
recent work of the author and collaborators P, [2] and [3]. Under the RWA (Rotating Wave 
Approximation) and some resonance conditions the problem is reduced to the evaluation of the 
exponential (i.e. finite time evolution operator) e~ ltH with the Hamiltonian 



H = 
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h 2 3 


h 2 4 











hu 




\ h u 


h 2i 


h 3 4 





) 



where hij are real coupling constants between the atom and laser fields. Here we have changed 
some notation from the previous one in [2] (like <-> for convenience. A generic 

algorithm to calculate e~ ltH based on eigenvalues of the Hamiltonian H has been given by us 
in [2], although the actual execution is not so simple. 

In this paper we revisit the problem from a different point of view. Let us introduce 
Makhlin's theorem. The isomorphism 

SU{2) ® SU(2) ^ SO(4) (^ su{2) ® 1 2 + 1 2 ® su{2) S so(4)) 

is one of the well-known theorems in elementary representation theory and is characteristic of 
four dimensional Euclidean space. In [1] Makhlin gave it the adjoint expression explicitly: 

F : SU{2) <g> SU{2) — ► SO (4), F{A ® B) = Q\A ® B)Q, 

with some unitary matrix Q G U(4). As far as we know this is the first time that the map was 
realized by the adjoint action. See also [5], where a slightly different matrix R has been used 
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in place of Q. This R (in our notation) is interesting enough and is called the magic matrix by 
Makhlin, see also [Sj and [7j. 

As the Hamiltonian H above is real symmetric (e su(4)) (not anti-symmetric (e so(4))) 
it is not obvious whether the magic matrix R could be applied to it or not. As we will show 
presently, a fairly wide sub-class of such Hamiltonians can be explicitly exponentiated with 
the help of the magic matrix and an additional similarity transformation. For the generic 
Hamiltonian H, we provide an approximate result of explicit exponentiation. 

Here we note that there is some overlap between our work and [7j. However, the methods 
given in [7] are quite varied and rather complicated. On the other hand ours rely on a consistent 
use of the magic matrix, so we believe our method could be easily digested by general readers, 
though the scope might be limited. 

As a closely related subject, we discuss the Baker-Campbell-Hausdorff (B-C-H) formula. It 
is one of the fundamental theorems in elementary Linear Algebra (or Lie group): 

e A e B = e BCH(A,B) . B CH(A, B) = A + B + 1 [A B] + J_ { p, B], B] + [A, [A, B]]} + • ■ • , 

where A and B are elements of some algebra. See for example the textbooks [H], or [TU]. A 
closed expression for the B-C-H formula for SU(2) is quite well-known. Similar closed expression 
is obtained for 5*0(4) (namely, for A,B e so(4)) by making use of the magic matrix, see [2]. 

Here we address the problem of explicit summation of the right hand side of the B-C-H 
formula for certain lower dimensional matrices, in particular, those related to su(4). 

In this paper we treat matrices of type 



H = 



( 


hi2 





h 14 









h 2 3 










h 2 3 





h 34 




\ hu 





h 3 4 





J 



which are very important in quantum optics, quantum computation, etc. (see [T]). We will 
derive its exact exponential form e~ ltH and will present a closed form expression of the B-C-H 
formula for two of them. 

2 Magic Matrix 

In this section we introduce appropriate concepts and notation together with a brief review of 
the results in [5] within our necessity, which look slightly different from the original work of 
Makhlin in @]. 
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The 1-qubit space is C 2 = Vectc{|0), |1)} where 
Let {o"i, o"2, cr 3 } be the Pauli matrices acting on the space 

*-(:;)■ «)■ — (o -0 - i2= (j 1 1- 

Next let us consider the 2-qubit space. Now we use the tensor product notation which is 
different from the usual one. That is, 

C 2 ®C 2 = {a ® b | a, 6 G C 2 }, C 2 ®C 2 = A iQj ® | a j} bj G C 2 , A,- G C, A; G Nj = C 4 . 
Then the 2-qubit space is 

C 2 §C 2 = Vect c {|00), |01>, |10), |11)}, 

where \ab) = \a) <g> |6) (a, b G {0, 1}). The Bell basis {|*i), |* 2 ), 1*3), 1*4)} defined by 

|*0 = -L(|00) + |11)), |* 2 ) = -L(|01) + |10)), |* 3 ) = ^(|01)-|10», |* 4 ) = -^(|00)-|11)) (3) 

plays an important role in various context. 

By H (2; C) we denote the set of all traceless hermitian matrices in M(2; C). It is well- 
known that they are spanned by the Pauli matrices 

H (2; C) = {a = a\0\ + a 2 a 2 + a 3 a 3 \ a±, a 2 , a 3 G R} 

and H (2; C) ^ su{2) where su{2) = £(31/(2)) is the Lie algebra of the group 5*7(2). 

The isomorphism is realized as the adjoint action (the Makhlin's theorem) as follows 

F : SU{2) <g> SU{2) — > SO{A), F(A <g> B) = R\A ® B)FL 

where 



R = (|*i), -i|* 2 >, -|* 3 ), -i|* 4 )) = i 



(4) 



/ 1 -i \ 

— i — 1 

-i 1 

\ 1 i J 

Note that the unitary matrix R is a bit different from Q in [1]. 

Let us consider this isomorphism at the Lie algebra level because it is in general easier than 
at the Lie group level: 
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£{SU{2) ® SU(2)) 



./ 



exp 



SU(2) <g> 5?7(2) 



£(,50(4)) 



exp 



50(4) 



Since the Lie algebra of SU(2) ® 37/(2) is 

£(5t/(2) <g> ST/(2)) = su(2) <g> 1 2 + 1 2 ® su{2) = {i(a ®l 2 + l 2 ®b)|a, be H (2; C)} 
we have only to examine 

f(i(a ® 1 2 + 1 2 ® b)) = 2i? f (a <g> 1 2 + 1 2 ® b)R E £(50(4)) = so(4). 
If we set a = Y^=i OjCTj and b = Y^=i bj a j then the right hand side of © reads 



«i2 f (a® 1 2 + 1 2 ® b)i? 



Conversely, if 



,4 





-(a x + 6i) 
-(a 2 - b 2 ) 



( f 12 

~fl2 
~ /l3 — /23 



ai + 6i 


-(a 3 - 6 3 ) 
a 2 + 6 2 

/l3 /l4 \ 
J23 hi 
/34 



a 2 - 6 2 
a 3 - 63 


-(ai - 61) 



A3 + &3 \ 

-(a 2 + 6 2 ) 
ai — 61 




e so(4) 



then we obtain 



V -/l4 "/24 "/34 7 

i?A_R t = i(a ® 1 2 + I2 ® b) 



with 



/l2 + /34 . /l3 — /24 . /l4 + ^23 

a = ai(7i + a 2 a 2 + a 3 cr 3 = <j\ H a 2 H 



^3, 



UI ,7 ,7 f 12 — fu fl3 + /24 . /l4 — /23 

b = b\G\ + 6 2 <T 2 + 63CT3 = <j\ a 2 H a 3 . 



(5) 



(6) 



(7) 



(9) 



It is very interesting to note that a and b are the self-dual and anti-self-dual part of the matrix 
A, respectively, under the Hodge ^-operation defined by {*F)ij = ~ Yl\i=i e ijkiFki- Here e^i is 
the totally anti-symmetric tensor with ei 2 34 = 1. 
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3 The Exponential of Matrices in Four Level System 



This section provides the main results of the paper, explicit calculation of the exponential of 
certain matrices. The generic Hamiltonian treated in this paper is a real symmetric 4x4 matrix 
with vanishing diagonal elements 



H 



( 


h-12 


hi3 


h 14 









h 2 3 


h 2 4 




hi3 


h 23 





h 3 4 




\ h 14 


h 2 4 


fl34 





) 



hij G R,. 



(10) 



The general method is quite simple: evaluate the evolution operator (matrix) U(t) = e ltH by 
unitary conjugation of the Hamiltonian H in terms of the magic matrix R: 



U(t) = e~ itH = RR)e~ itH RR) = Re~ itB?HR R). 



The conjugated Hamiltonian R^HR reads explicitly 



R^HR 



1 



/ 2h 1A 

i(h 12 + hu + h 2 4 + h 3i ) 

-hi2 + /li3 - /l 2 4 + ^34 

o 



-i{h\2 + hi 3 + h 2 4 + h 34 ) -hu + h 13 - /i 24 + h 34 

2h 23 

-2^23 

hi2 + h 13 - h 2 4 - h 34 -i(hi2 - h 13 - h 2 4 + h 3i ) 



\ 

h\2 + h\ 3 — fl24 — h 3 4 

i{h\2 - h 13 - h 2i + h 3i ) 
-2hu 



I 

(12) 



whose structure is better displayed in the following tensor product notation 



R ] HR 



'13 



4 24 



-CTi 



/ii4 + h 



23 



-a 3 



h + h 



his + h 24 hu — h 2 3 
cr 2 H C7 3 



. h\2 + ^34 _ 
H ~ 0"3 ® &2 



h 



12 



h 



31 



-<J\ QS> cr 3 . 



(13) 



2 ° 2 

Obviously the first two terms commute with each other and the last two terms U3 ® a 2 and 
<J\ ® 03 also commute. The right hand side of (TIB"]) is a summation of two blocks not commuting 
each other. 



A few comments are in order. 

(1) In the calculation above we have used R)HR in stead of RHR) in (JTj). This point is 
important (we leave it to the readers to contemplate why it is so). 

(2) The space of the entangled states in the two-qubit system is SU "(4) / ' SU ~{2) ® 577(2), which 
is identified with the homogeneous space SU(A) / SO(4) because of the isomorphism SU{2) <g> 
577(2) = 50(4). Then it is well-known that 



[A e SU(A) I A T = A] = {A T A | A e SU(4)} = SU{A)/SO{A), 
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K 



diag(7ii, h 2 , h 3 ,h 4 ) + H, hi + h 2 + h 3 + h A = 0, (14) 



see for example [5]. By use of the expression SU(4) 3 A = e lK with K G -Ho (4; C) we have 

SU{A)/SO(A) = {e iK | K e H (4; R)} , 

where H Q (A; R) is the set of all traceless symmetric matrices in M(4; R). It is of dimension 9 
(see flUD), as expected: 15(sw(4)) - 6(so(4)) = 9. 
Since this K is written as 

/ h h 12 h 13 h 14 \ 
hi 2 h 2 h 23 h 2A 
hi 3 h 23 h 3 /i 34 
\ h\4 h 2 4 h 3A h A ) 

the evaluation of the evolution operator U (t) = e~ ztH is deeply related with the study of the 
entangled states in the two-qubit systems. 

For later convenience the conjugated matrix of K is given here: 

Dtr , D / h X3 -h n h u + h 23 \ / h 13 + h 24 h u - h 23 
WKR = I a x + a 3 1 <8> 1 2 + 1 2 <g> I a 2 + a 3 

, h\ 2 + /l 3 4 /Zi2 — /i34 ^ 

H cr-3 ®cr 2 - cti <8> a 3 

h - h - h 2 + h 3 h + h!-h 2 -h 3 h - h\ + h 2 - h 3 
H cr 3 <g> cr 3 H (Ti <g) cr 2 H cr 2 ® 0"i. 



(15) 



3.1 Approximate Result 



(17) 



We do not know yet how to exponentiate the generic form of the (i?-conjugated) Hamiltonian 
ffTHj) . Our first main result is an approximate one for the generic case: 

U{t)n 11^)11^)11^)11 A {t), (16) 

where each factor in ffTHl) can be exponentiated exactly, 

U^t) = R {exp (-it (^J^a, + ^^a 3 )) ® l a } rf, 
f/ 2 (t) = i2 {l 2 ® exp (-if (^+J^ 2 + ^±^*<r 3 )) } (18) 
Z7 s (t) = i?exp (-it ^ 12 ^ 3 V 3 ® ^ Ht, (19) 

*7 4 (t) = Rexp (-it ( y - hl2 ~ 2 hM a 1 ® a 3 Yj Rl (20) 

As remarked earlier, U\ and £/ 2 commute with each other and U 3 and U4 also commute but the 
other pairs do not commute. Hopefully the approximation is not so bad. The non-exactness 
arises from the non-commutativity. 
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3.2 Special Exact Result I 

If h\2 = h 3i = in ffT3"]) the troublesome non-commutativity disappears. Then we have the 
simple form 

Otero f h n-h 2 4 h u + h 23 \ ( h 13 + h 24 h 14 - h 23 \ 

WHR = I cri H a 3 J <g) 1 2 + 1 2 ® I o 2 H a 3 J (21) 

with 



/ 





^13 


ft<14 










^23 


^24 






^23 










\ ^14 


^24 








) 



H 



Then we have the exact form 

U{t) = U 1 {t)U 2 {t) = U 2 {t)U 1 {t) 



{22) 



(23) 



with Ux{t) in (TUD and t/ 2 (t) in (fig]). 

Though this Hamiltonian is of very special form there is some application shown in the 
following section. 



3.3 Special Exact Result II 

The Hamiltonian that we really want to study is 



H 



( 


h.12 







\ 


h 2 





^23 










h 23 





^34 




V h u 





h 3 4 





I 



(24) 



which is more general than the one discussed in recent papers by the author [T] [2] and [3], 
which is obtained by setting h\ 4 = 0. This restricted case was also discussed in [7]. At first 
sight this Hamiltonian looks rather different from ( J22|) . which can be exactly exponentiated. 
However, it can be reduced to the form of f)22p by a similarity transformation in terms of the 
exchange (swap) matrix S: 



SHS 



( 





hi2 


h 14 










h 23 


hu 




hi2 


h 23 










\ h u 


h 3i 








) 



(25) 
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with 



S = ~(h ® h + O"! <8> O"! + 0"2 ® 0"2 + 0"3 ® 0"3y 



/ 1 \ 

10 

10 

\ 1 / 



5 = s T = s- 1 



Therefore we have 



hi? — ho, 



hi * + h-) 



^^a 1 + ^i^a 3 )®l 2 + l 2 ®(^^^ (26) 



U^t) = SR 



|exp ^— it f 



hl2 — ^34 . hi4 + h 2 3 

^ 0i H ;r 



2 2 V V 2 2 

The evolution operator 

C/(t) = e"^ = SRR^Se~ itH SRR^S = SR e - itR " SHSR R^ S 

takes an exact factorized form 

[/(*) = U x {t)U 2 {t) = U 2 {t)Ui{t) (27) 

where 

r / / f, f, i /, \\ 

(28) 
(29) 

Obviously the two factors commute and they belong to two independent 577(2). Compare this 
result with that of [7j. 

4 B-C-H Formula for a class of matrices in SU(4) 

In this section we give a closed expression to the B-C-H formula for the set of unitary matrices 
{e zH | if is a type of fpH|) } by use of the results in the preceding two sections. Let us prepare 
some notation for simplicity. 

For two (real) symmetric matrices A, B 



U 2 (t) =SR{1 2 ® exp ( -it [^±^a 2 + ^^a 3 ) 



.4 



/ o h o u \ 

fi f 2 

o h o h 

\ U o h o / 



B 



( g x g 4 \ 
9i g 2 

\34 53 0/ 



(30) 



we can set 



R*(SAS)R = ai ® 1 2 + U ® a 2 , R\SBS)R = b x ® 1 2 + 1 2 ® b 2 
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Let us represent these four su(2) elements in terms of the R 3 vectors as in (135]) : 



ai 
bi 



fl — fl . /4 + J2 

— 2 — ff i + — ^ — CT3 ' a2 

01-03 . 04 + 02 , 

— 7, — °i H ^ — °"3, b 2 



fl + J3 . J 4 — 

— &2 H 7> 0"3> 



2 

01+03 



-a 2 



2 

04 - 02 



-CT 3 



with 





/ /1-/3 
2 






/ 91-33 
2 






f 






f 


\ 


ai = 







, bi = 







; a 2 = 


/1+/3 
2 




, b 2 = 


31+33 
2 






1 fr+fu 
\ 2 






L 94+32 
\ 2 


J 




l /4-/2 
\ 2 


J 




, 34-32 
\ 2 


/ 



Obviously a x and bi belong to the same su(2), whereas a 2 and b 2 belong to the other su(2). 
Next we introduce a pair of three real parameters (ai,/3i,7i) and (a 2 ,/3 2 ,7 2 ) defined for each 
pair of R 3 vectors (a"i, bi) and (a 2 , b 2 ): 

ai=a(ai,bi), p 1 =p(a 1 ,b 1 ), 71 = 7(^1, bi), 
a 2 = a(a 2 ,b 2 ), (3 2 = /3(a 2 , b 2 ), 72 = 7(^2, b 2 ). 

See ( 1371) in the Appendix or the paper [TT] for the definition of these parameters. By combining 
the exact B-C-H formula for su(2) for each pair (ai,bi) and (a 2 ,b 2 ), we obtain 

e iA e iB = SRR^Se lA SRR^Se iB SRR^S 
= SRe iRtSASR e iRtSBSR R ji S 

= SR (e iai ® e ia2 ) (e 4bl ® e iba ) R^S 
= SR (e iai e 4bl ) ® (e ia2 e ib2 ) R)S 

= 5'J^e i ( aiai+ ^ lbl+ '>' 1 f [ ai ' bl l) (g) e *( Q2a2 +/3 2b2 +72§ [a2,b2])^f c> 

= ^^e i ^ aiai+ ^ lbl+7l ^ ai ' bl ^® l2+l2 ®^ a2a2+ ^ 2b2+72 ^ a2 ' b2 ^J R) S 



iSR{ (aiai+ J 3ibi+7i|[ai,bi])®l2+l2®(Q2a2+/32b 2 +72f [a 2 ,b 2 ])}fltS 



_ e iBCH(A,B) 



The desired closed form of the B-C-H formula reads 
BCH(A,B) 

3> I2 + I2 < 



(31) 



ST? I ( aiai + /?ibi + 71 - [ai, bi 

/ _0_ (12) (13) (14) \ 

(12) (23) (24) 

(13) (23) (34) 
V (14) (24) (34) J 



ct 2 &2 + /?2b 2 + 72^ [a 2 , b 2 ]^J J R*S 



(32) 
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whose entries are 

n \ fi ~ fa . a 9i — 93 . fl + fa . a 9l + 93 

(12) = at — - — + ft — - — + a 2 — - — + ft — ^ — ' 

/, Q x ■/ ( h- h9i + 92 U + f2 9l-93\ ffl + f3 94 ~92 fi~f2 91 + 93 

(13) = i <^ 7i — — - 72 



f-iA\ fa + fa , a 94 + 92 fa~ fa . a 94 — 92 

(14) = on — - — + ft — - — + a 2 — - — + ft — - — , 

/ OQ \ fa + h , a 94 + 92 fi — $2 n 94 — 92 

(23) = a x — - — + ft — a 2 — ft — - — . 

foA \ ■/ / h ~ h 94 + 92 fa + f2 9i~93\ . (h + h94 - 92 fa-f2 9i + g3 

(24) = i i 7l — — +72 



2 2 2 2 J ' \ 2 

/ Q/ ,N fl~f3 a 9l — 93 . fl + fa . Q 9l + 93 

(34) = —ol\ — ft — h a 2 — h ft — - — . 



Here (12) is the complex conjugate of (12). 

This is another main result of the present paper. 



5 Discussion 

In this letter we addressed the problem of explicit exponentiation of the generic Hamiltonian in 
the general four level system. With the aid of the magic matrix, we obtained an approximate 
result for the generic case and some exact results for certain restricted forms of the Hamiltonian. 
It is a good challenge to derive the explicit form of the evolution operator for the most generic 
Hamiltonian in the four level system. 
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Appendix: B-C-H Formula for SU(2) 

In this Appendix we recapitulate the closed form expression of the B-C-H formula for 577(2) 
in the 2x2 representation, see [H] for more details. 

First of all let us recall the well-known exponentiation formula: 

e i(xa 1+ ya 2+ za 3 ) = ^ ^ + + + Z(J \ (33) 

r 



i— {x + iy) cosr — i^-z 



(34) 
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where r = yj x 2 + y 2 + z 2 . This is a simple exercise. 

For the group SU(2) it is easy to sum up all terms in the B-C-H expansion by using the 
above exact exponentiation formula. Because of the obvious relation su(2) = so(3), three 
dimensional vector notation is quite useful for su(2). For 2x2 hermitian matrices, 

X = xiai + x 2 a 2 + x 3 a 3 , y = V\(J\ + y 2 o 2 + Vzoj, e H (2, C) 

we associate R 3 vectors. Namely, we set 

fx, \ 



X 



x 2 

\X3 ) 



Y 



(yi\ 

\f3 / 



and 



x ■ y = xxyi +x 2 y 2 + x 3 y 3 = Tr(XF)/2, |x| = a/xoc, |y| = y/y-y- 
Then the commutator of X and Y corresponds to the vector product x x y: 



/ 



x x y 



x 2 y 3 - x 3 y 2 
x 3 yi - xiy z 
\ x\y 2 - x 2 yi ) 



(35) 



Now we are in a position to state the B-C-H formula for SU(2): 



e iX e iy = e iZ : Z = aX + (3Y + j-[X, Y], 



where the real coefficients a, (3 and 7 are defined by 

= «l x > yj = : -n — — ' ^ = ^( x > y) 

p m 

«in~ n «in IyI «in Ivl 

7 = 7(x,y) 

with 



sin 1 p sin 


X 


COS 


|y| 


P 

sin" 1 p sin 


X 


X 

sin 


y| 


P 


X 


|y 





sin 1 p cos 


X 


sin 


y| 


P 


y 





P 2 = p(x,y) 



sin 2 |x| cos 2 |y| + sin 2 |y| 



sin 2 


X 


sin 2 


y| 




X 


2 


y 


2 



x-y) + 



2 sin 


X 


cos 


X 


sin 


y 


cos 


y| 




X 


|y 





(x-yj 



(36) 



(37) 



The proof is not difficult, so it is left to the readers. 
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